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Method of Normalizing Color Information, Method of Exchanging Color 
Information using the Same Method and Device Using the Same Methods 

5 Field of the Invention 

The present invention relates to a method of normalizing color 
informatioruand a method of exchanging color information as well as a device** * ^ 
implementing this method. Those methods and the device can realize device- 
independent color reproduction in a color management system— managing color 
10 information on a network to which input/output image devices such as a camera, 
scanner, monitor, printer and the like are hooked up. 

Background of the Invention 
Advancement of the digitalization in image devices and development of 
15 network technology particularly in the Internet promote to ^prtrlrmTrr q cross- 
media-system. This system allows arbitrary exchanging of image information 
by coupling various image devices(puch as a scanner, printer and displas^eaclT 
other on an open system. Businesses such as E-commerce and Telemedicien 
medical service requiring strict color reproduction have taken off thanks to this 
20 system. 

The open system uses various devices in various viewing environments. 
* ]Sffi£feaaiaflie and hardware, which are independent of both the device colorimetoric 
characterization and viewing environment^allow color reproduction of a 
certain quaht^aam^v ^aa^m^n^^^ n rxri^ m - 

25 A camera or a scanner (a source device) should deliver the color 

information it takes in to the open system exactly. On the other hand, a display 
or a printer (an output device), which outputs the color information it receives, 
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should render the information exactly. 

For instance, when a camera takes color information exactly, <smA if a 



display caasaaitl processrfthe information inadequately, then the exact color 
information cannot be displayed. An accuracy of color reproduction by this 
5 entire system is thus lowered. 

For instance ^4fa»§KSDgih E-commerce vi^Internet, a customer views goods 
on a display in different colors from actual ones. 

However, when a color displayed area (gamut) taken by a source device is 
q larger than the gamut of an output device, the source device should change the 

m 10 color information to fit within the gamut of the output device. This 
.7j manipulation is referred to as *?kmut mapping". The color outside the color 

gamut of the output device cannot be reproduced; however, when a color image 

• , n ■ ftttlfclfcH 

w is reproduced, an image painted by colors as^close as the original image is 

I- required. Therefore, gamut mapping should be optimized, so that a reproduced 

!=& 15 image appears as ffi&to&aaiH^ as possible. 

ill 

O In the color management system (CMS), the following elements should be 

,5ST optimized because color information must be exchanged correctly between 

different image devices as well as viewing environments^ 

(1) Input/output characteristics (relation between values of device- 
20 driving-signal and colorimetric values) of respective image devices is 
standardize^ 

The principle of the CMS is to establish input/output characteristics, i.e^a 
relation between a video signal (e.g. RGB signals) driving an image device and 
colorimetric value (e.g. CIEXYZ of XYZ display system or CIELAB of L,a,b, 
25 display systenLboth are authorized by Commission International de I 1 Eclairage 
(CIE)). Then the colors from the source device are converted into device- 
independent color data so that any device^ can handle the color data through a 
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given procedure. Then first of all, the input/output characteristics of the 
input/output devices linked to an open system should be standardized. 

International Electrotechnical Commission (IEC) is modeling the 
input/output characteristics of image devices such as an LCD, PDP, scanner, 
5 digital camera^pwaafea&^iifcd drafting a standard specification. 

(2) Input/output characteristics of viewing-environment-dependent 
visual system^s^tandardized. 

Tri-stimulus value of CIEXYZ is a psychophysical quantity produced by 
multiplying a physical quantity of light and a spectral-response characteristic of 

10 human eyes (when it is expressed with a wave-length function, it is referred to 
as color matching function). Two colors having the same XYZ thus appear the 
same. However, when the color matching function — expressing spectral- 
response characteristic of the human eyes — changes depending on viewing 
environments and human chromatic adaptation is in different viewing 

15 environments, the two colors having the same XYZ may appear differently. 
This phenomenon is referred to as, in general, "different color appearance". An 
objective of the CMS is to reproduce, in any viewing environments, the color the 
output device receives as a color appearance as same as the color information 
the source device sent out. 

20 It is thus necessary to standardize the influence of viewing environment 

which the input/output characteristics of a human visual system is subjected to 
in order to find a counterpart color (corresponding color) of the color XYZ which 
appears the same under a different lighting condition on a source device and an 
output device. 

25 The CIE specified a color appearance model CIECAM97s in TCI-34 

(Testing Color Appearance Models, Report on CIECAM97s, April 1998) where 
the model considers the brightness of background and the XYZ values of a white 
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point influencing the human chromatic adaptation. 
(3) Gamut mapping is optimized. 
A method is studied to send out color information by compressing the 
entire gamut of a source device in order to fit the gamut within a smaller gamut 
5 of the output device. 

Many engineers have proposed various methods of gamut compression, 
and the CIE is now trying to standardize the method. 

Regarding the method of gamut compression, GCUSP developed by Jan 
Morovic has been well known; however, this method has the following problems: 
10 The GCUSP is a gamut compression method to minimize the lowering of chroma, 
however, it does not consider the three attributes of color reproduction 
comprehensively, i.e. lightness, chroma and hue after performing the gamut 
mapping. The color reproduction does not solely depend on chroma but the 
lightness and hue are also the critical factors. 
15 The method of gamut compression cannot reproduce the same color as a 

color appearance or colorimetric value when a source device has different gamut 
from an output device. As a result, a problem occurs that an agreement 
between colorimetric values or color appearances cannot be applied to as a 
design target. 

20 A method of agreeing the color appearances is studied. One of effective 

methods is to control colors based on an agreement of color names. 

When the gamut mapping is performed, it is impossible to render the 
same colors on respective devices, thus it will be a target to exchange color 
information between respective devices in order to show the difference in overall 
25 image as little as possible. 

To achieve this target, a concept of categorical color perception is 
introduced. A human being has an ability of perceive a color and another 



ability of grouping some colors into a color category. For instance, light green', 
? dark green 1 and ! bluish green' appear to be different when they are placed 
sequentially; however they are grouped into one category 'green 1 . 

If a red flower on a CRT display is printed ! a red flower 1 in the same 
category on a sheet of paper, viewers do not mind subtle differences and get an 
impression that the two images are the same; however, when it is printed in a 
different color category such as ! a yellow flower' or f a flower of orange color', the 
viewers get an impression that the two images are different. 

Further, when a flower's red has more intensity than an apple's red on a 
CRT, it is preferable to see the flower's red with more intensity than the apple's 
red on a printed matter, and it is not preferable to see the flower's red with less 
intensity than the apple's red. 

It is effective to exchange color information thereby agreeing color names 
between a source device and an output device, and yet, preserving a perceptive 
quantity relation of a color such as chroma and the like in order to minimize a 
difference of entire images in gamut mapping. 

Let's consider colors having a color difference between two points existing 
in a color space of a source device. When the color difference is strictly 
preserved and the two colors are mapped in an output device, the difference of 
the two colors from the viewpoint of the color difference is the same on the 
source device and the output device. However, when a color categorical control 
is not built in the CMS, the two colors in the same color category on the source 
device may be mapped at the spots striding over a boundary between color 
categories and reproduced in two colors having different color names. 

A method of color adjustment based on an agreement of color names is 
disclosed in the Japanese Patent Application Non-examined Publication No. 
H08-191400. This method is manipulated by an operator when he or she 
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adjusts a color rendered on a display so that the color can agree with a given 
name of a target color. 

To be more specific, the colorimetric value of each pixel of an image is 
expressed in lightness (L), chroma (C) and hue (H) in a polar coordinate system. 
Then the differences of each element between a reference color and a target color, 
i.e. AL, AC and AH are added to the color before adjusted, and the color after 
adjusted is obtained. Among the colors after adjusted, the color having color 
names different from those of target colors are re-adjusted with respect to the 
differences AL, AC and AH. In re-adjusting, AL/n, AC/n and AH/n are 
used (n: numeral around 10). If the re-adjustment cannot find the agreement 
with the color name of the target color, change the 'n' and repeat the adjustment 
of the differences until the color name finds the agreement. 

However, the color adjustment method discussed above has still have the 
following problems: 

First, differences ofAL, AC and AH are changed repeatedly based on 
trial and error base until the color name find the agreement. This procedure 
incurs increasing of work time, and yet, an image after the adjustment tends to 
be dispersed. 

Second, the names of colors adjusted with the differences of A L, AC and 
AH differ from those of the target colors, and they are re-adjusted with the 
differences of AL/n, AC/n and AH/n. In this case 'n' is determined by trial 
and error process. The trial and error process incurs increasing of work time, 
and the image after the adjustment is vulnerable to dispersion. Further, when 
a plurality of adjusted colors are available, determining 'n' by trial and error 
process may reverse the positions of the adjusted colors in color space, or 
produce a color-skip. These phenomena are viewed as a pseudo outline and 
degrade the color reproducibility. 
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Summary of the Invention 
The present invention addresses the problems discussed above, and 
provides a method of exchanging color information and a device using the 
5 method. The method is to group the CIELAB color spaces both of a source 
device and an output device into a plurality of color categories, and produce 
normalized color space by using gravity vectors calculated for each color 
category of respective devices as well as variance covariance matrix expressing 
Q extension from a point of gravity as a center,. Then control points of the same 

5 10 color categories provided in the respective normalized color spaces are mapped, 
*j* and an agreement of color categories between the two normalized color spaces is 

found, then color perceptive quantities such as chorma and the like of the two 
corresponding colors in the normalized color spaces are preserved, whereby the 
color information is exchanged. This method produces difference as little as 
j=£ 15 possible in overall impressions of the images of respective devices, 
p A color-information-exchanging-device of the present invention, in an 

w apparatus transmitting the color information from a source device to an output 

device, comprises the following five elements: 

(1) a color-naming-database for the source device, storing: 
20 a colorimetric value cox*responding to a color name a viewer 

responds to the color stimuli; 

a gravity vector of respective colors calculated from the colorimetric 

values; 

a variance covariance matrix expressing an extension of a color 
25 category calculated from the colorimetric values, 
by rendering color stimuli on the source device, 

(2) a color-naming-database for the output device, storing: 
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a colorimetric value corresponding to a color name a viewer 
responds to the color stimuli; a gravity vector of respective colors calculated 
from the colorimetric values; a variance covariance matrix expressing an 
extension of a color category calculated from the colorimetric values; 
5 by rendering color stimuli on the output device; 

(3) a normalizing section of color information for the source device for 
calculating a normalized distance between the gravity vector and the 
colorimetric value of the source device based on the variance covariance matrix; 

(4) a mapping control section for converting the normalized distance to 
10 target normalized distance which is a normalized distance in a color space of the 

output device with a coefficient by comparing the gravity vector read out from 
respective databases with covariance matrixes thereby finding a conversion 
coefficient from the comparison difference; and 

(5) a target normalized distance searching section for searching the 
15 color space of the output device for a colorimetric value at the target normalized 

distance from the gravity of respective colors. 

This color-information-exchanging-device avoids the difference in overall 
impression between a reproduced image on the output device and an original 
image on the source device even when the gamut of the source device differs 
20 greatly from that of the output device. 

Brief Description of the Drawings 
Fig. 1 is a block diagram of a color management system incorporating a 
first exemplary embodiment of a color-information-exchanging-device of the 
25 present invention. 

Fig. 2(a) shows a gamut of a CRT display and a printer. 

Fig. 2(b) shows another gamut of a CRT display and a printer. 



Fig. 3(a) shows a color distribution of a red category and a pink category 
displayed on a CRT display. 

Fig. 3(b) shows a color distribution of a red category and a pink category 
printed by a printer. 

Fig. 4 shows the color region indicated by Kelly. 

Fig. 5 shows the chromaticity diagram made by Mac Adam 

Fig. 6 illustrates a relation between intersection lines of Mac Adam's 
chromaticity diagram and boundary lines of Kelly's color region. 

Fig. 7 illustrates an inside structure of color-information-exchanging- 
device 101 in accordance with the first exemplary embodiment of the present 
invention. 

Fig. 8 illustrates an inside structure of gravity mapping controller 704 in 
accordance with the first exemplary embodiment of the present invention. 

Fig. 9 illustrates an inside structure of a dynamic range mapping 
controller 705 in accordance with the first exemplary embodiment of the present 
invention. 

Fig. 10 illustrates a system used in an experiment of color naming. 

Fig. 11 shows a procedure of mapping into printer's space with preserving 
color categorical characteristics of the color space of a CRT display. 

Fig. 12 shows a procedure of searching the color space of an output device 
for a colorimetric value having a target normalized distance. 

Fig. 13(a) shows a cusp of gamut of a source device and that of an output 
device. 

Fig. 13(b) shows gamut boundaries of the source device and output device, 
the boundaries are indicated by a plane of chroma C — lightness L. 

Fig. 14 is a block diagram of a color-information-exchanging- device in 
accordance with a second exemplary embodiment. 
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Fig. 15 illustrates processes of a main mapping and a pre-mapping when 
lightness mapping from a CRT display to a printer is performed in accordance 
with the second exemplary embodiment of the present invention. 

Fig. 16 shows cross sections of CRT display's gamut and printer's gamut. 
5 Fig. 17(a) is a projection drawing of respective cusps of a CRT display and 

a printer of CIELAB space projected to a plane of a — b. 

Fig. 17(b) is a projection drawing of respective cusps of a CRT display and 
a printer of CIELAB space converted to a polar coordinate projected to a plane of 
n L_H - 

'P* 10 Fig. 18(a) illustrates a case where the lightness L s of a source device is 

9 smaller than cusp's lightness L s cusp mapped mapped on gamut surface of an output 

=C device. 

Fig. 18(b) illustrates a case where the lightness L s of the source device is 
greater than cusp's lightness L s cusPtmapped mapped on gamut surface of the output 
15 device. 

Fig. 19(a) shows a distribution of color category on the gamut surface of a 
□ CRT display. 

Fig. 19(b) shows a distribution of color category on the gamut surface of a 
printer. 

20 Fig. 20 shows mapping result of lightness of CRT display's cusp. 

Fig. 21 illustrates an inside structure of a pre-mapping section in 
accordance with the second exemplary embodiment of the present invention. 

Fig. 22 illustrates an inside structure of a main mapping section in 
accordance with the second exemplary embodiment of the present invention. 
25 Fig. 23 illustrates an inside structure of a gravity control point mapping 

section in accordance with the second exemplary embodiment of the present 
invention. 
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Fig. 24 illustrates an inside structure of a surface control point mapping 
section in accordance with the second exemplary embodiment of the present 
invention. 

Fig. 25 illustrate chroma C Ci S) maxd and chroma C di maXj t at finding gamut 
shape control coefficient w c . 

Fig. 26 is a block diagram of a color management system introducing a 
color information exchanging device in accordance with a third exemplary 
embodiment of the present invention. 

Description of the Exemplary Embodiments 
Exemplary Embodiment 1 

The first exemplary embodiment of the present invention is demonstrated 
hereinafter with reference to the accompanying drawings. 

Fig. 1 is a block diagram of a color management system (CMS) 
incorporating the first exemplary embodiment of a color-information- 
exchanging-device of the present invention. 

For simplifying the description, an image on a CRT display is printed out 
by a printer in Fig. 1. 

The CRT display is used as a source device and the printer is used as an 
output device. In Fig. 1, color-information-exchanging-device 101 converts the 
color information (colorimetric value XI, Yl, Zl) displayed on CRT display 102 
into corresponding color information (colorimetric value X2, Y2, Z2) of printer 
104. Device 101 then outputs the converted information, where the 
corresponding color information appears the same as the color information of 
the CRT display. CRT display driving unit 103 supplies a drive signal (R,G,B 
signals) colorimetrically adequate into CRT display 102 based on display's 
input/output characteristics described in a CRT display profile. Driving unit 
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103 also converts the drive signal into tri-stimulus values (XI, Yl, Zl) — color 
information independent of input/output characteristics of display 102, then 
outputs the converted values to color-information-exchanging-device 101. 
Printer driving unit 105 converts tri-stimulus values (X2, Y2, Z2) 
colorimetrically adequate into a drive signal (C, M, Y, K signals) based on the 
printer's input/output characteristics described in printer's profile 107. Then 
driving unit 105 outputs the converted signal to printer 104. 

Before the color-information-exchanging-device of the present invention 
is demonstrated, some premises are described hereinafter. 

First, each gamut of the CRT display and the printer is described. 

Fig. 2 shows an example of each gamut of the CRT display and the 
printer. 

Fig. 2(a) is a projection drawing where the gamut is projected to the plane 
of a~ b of CIELAB's color space. 

In Fig. 2(a), inside of gamut boundary 201 is a gamut of the CRT display. 
Inside of gamut boundary 202 is a gamut of the printer. 

Fig. 2(b) shows cross sections of each gamut being cutoff at hue-angle 0 
degree in a cylindrical coordinate converted from the coordinate of CIELAB's 
color space. In the cylindrical coordinate, the horizontal-axis represents 
chroma (C) and the vertical-axis represents lightness (L). In Fig. 2(b), inside of 
gamut boundary 203 is a gamut of the CRT display and inside of gamut 
boundary 204 is a gamut of the printer. Fig. 2(a) and Fig. 2(b) indicate that the 
CRT display has a greater gamut than the printer, and when an image on the 
CRT display is output to the printer, gamut mapping is needed. 

Second, a color temperature is described hereinafter. A typical white 
point's color temperature of the CRT display is 9300K. The color temperature 
of a paper white of a printer can be adjusted at the same temperature as that of 
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the CRT display at the printing site. However, a standard viewing condition 
ISO 13655 specifies the viewing temperature at 5000K. Therefore, the paper 
white of the printer shows a color temperature 5000K under the standard 
viewing condition. Actually, it is a rare case that a color temperature of a room 
5 lamp agrees with 5000K. 

CRT displays are often viewed under the room lamp, and the light of the 
room lamp is intermingled with self-emitting light from the CRT display, 
therefore, a gamut mapping dealing with an arbitrary color temperature is 
demanded. 

10 Regarding the color temperature, there are cases where the CRT display 

has the same color temperature as the printer, while the CRT display has a 
different color temperature from the printer. 

In general, an image displayed on the CRT display often does not match 
with an image printed by the printer. Therefore, a predictable mapping — the 
15 color rendered on a source device is predicted how the corresponding color 
appears on an output device — is required. 

Third, a color categorical characteristic is described hereinafter. 
Fig. 3 shows experimental data including color-category-names 
(categorical color naming) viewers give to the colors displayed on the CRT 
20 display or printed by the printer. 

In the CIELAB's color space, three-dimensional color vectors are 
produced at random, and the three-dimensional vectors within a reproducible 
range as test colors are shown to viewers. Viewers then select a name out of 11 
basic categorical color names, i.e. white, black, red, green, yellow, blue, purple, 
25 pink, orange, gray and brown. The experiment has been done in this way 
discussed above. 

In Fig. 3, selected red (represented by black points) and selected pink 
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(represented by white points) are plotted on the plane of a~ b of CIELAB's color 
space. Fig. 3(a) shows a result of the naming by displayed colors on the CRT 
display, and Fig. 3(b) shows a result of the naming by printed colors. 

When Fig. 3(a) and Fig. 3(b) are compared, there are differences in color 
5 categorical distribution even in a common gamut shared by both the devices. 
This means that when the bare CIELAB value is sent from the CRT display to 
the printer, a color-categorical mismatching occurs. 

In Fig. 3, color 301 displayed on the CRT display and color 302 printed 
have the same CIELAB values; however this color is perceived as 'pink 1 on the 
10 printed paper, and that color is perceived as 'red' on the CRT display. 
Accordingly, when a color displayed on a CRT display and a color having the 
same CIELAB value as the color displayed on the CRT appear on a sheet of 
paper by a printer and the printed colors have different color names. This 
phenomenon lowers color reproducibility. Thus the mapping should be 
15 controlled so that the two colors may have the same color names. This control 
is a key to improve the color reproducibility 

Fig. 4 shows the color region indicated by Kelly (reference document: 
Basic Color Engineering written by Mitsuo Ikeda, pllO) and Fig. 5 shows the 
chromaticity diagram made by Mac Adam (reference document: ditto pl37). 
20 The ovals are enlarged 10 times for easy understanding. 

This result leads to the following prediction: A group of colors, e.g. 
'green', 'yellowish green' and 'bluish green' which Kelly used, can be expressed in 
three-dimensional oval. 

Fig. 6 illustrates a relation between intersection lines of Mac Adam's 
25 chromaticity diagram and boundary lines of Kelly's color region. In Fig. 6, oval 
601 indicates a group of color which is perceived 'green'. Oval 602 indicates a 
group of 'yellowish green' and oval 603 indicates 'bluish green'. 
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Boundary line 604 indicates a boundary between 'green' and 'yellowish 
green', and boundary line .605 indicates a boundary between 'green' and 'bluish 
green'. Boundary line 606 indicates a limit of CRT display's displayable colors, 
and color region above the boundary are not displayable by the CRT display. 

Based on those premises, * the color-information-exchanging-device in 
accordance with the first embodiment is described hereinafter. In this case, the 
CRT display and the print have the same reference white, namely, the case 
where colorimetric mapping should be performed. 

Fig. 7 is a block diagram of the color-information-exchanging- device in 
accordance with the first exemplary embodiment of the present invention. 

In Fig. 7, color-information-exchanging-device 701 outputs colorimetric 
value of the color displayed on the CRT display as the colorimetric value of the 
CRT display. Color-naming-database 702 for the source device stores (a) a 
colorimetric value of the color name selected among the colors displayed on the 
CRT display by at least one viewer, (b) gravity vector of each color calculated 
from the colorimetric values, and (c) variance covariance matrix calculated from 
the colorimetric values and expressing an extension of a color category of each 
color. Color-information-normalizing-section 703 for the source device reads 
out colorimetric values of respective colors from database 702, calculates gravity 
vectors and variance covariance matrixes of respective colors, and stores the 
resultant calculations into database 702. At the same time normalizing section 

703 calculates a normalized distance between the CRT display's colorimetric 
value and the gravity vector based on the variance covariance matrix and 
gravity vector. The normalized distance undergoes gravity mapping controller 

704 and dynamic range mapping controller 705 and is converted to a target 
normalized distance, then fed into target-normalized-distance-searching section 
706. 
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Color-naming-database 707 for the output device stores (a) a colorimetric 
value of the color name selected among the colors printed by the printer by at 
least one viewer, (b) gravity vector of each color calculated from the colorimetric 
values, and (c) variance covariance matrix calculated from the colorimetric 



Target-normalized-distance-searching section 706 reads out colorimetric 
values of respective colors from database 707, calculates gravity vectors and 
variance covariance matrixes of respective colors, and stores the resultant 
calculations into database 707. At the same time, searching-section 706 
calculates a normalized distance between candidate colors for outputting and 
the gravity vector in the printer's color space, then finds a colorimetric value 
having a target normalized distance, and outputs the resultant values to color- 
information-output-section 708. 

Color-information-output-section 708 outputs the resultant colorimetric 
value as the printer's colorimetric value. 

Next, a calculation method to find a normalized distance of color- 
information-normalizing-section 703 for the source device is detailed. 

A color category can be shaped in an oval as discussed previously, thus a 
distance between the color categories is expressed with a normalized distance. 
The normalized distance Di is found through the following equation: 



values and expressing an extension of a color category of each color. 




(1) 



where X = color vector; 



Mi = gravity vector of color category T; 

Zi = variance covariance matrix of color category ! i'; 

Ei _1 = inverse matrix of the variance covariance matrix Si; and 
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Di = colorimetric distance from the gravity of the color category 'i' to color 
vector X, and a distance normalized by Si expressing an extension of a color data 
group of color category T. 

In CIELAB color space, X, fii and Si are denoted as follows. 

X=[ Xl x 2 x 3 f (2) 



(3) 



2r 



$11 $12 $13 

S 21 S 22 S 23 

■^31 ^32 ^33 
1 



— g - Mi)' ^ 2 < " M,)(* 2i - M 2 ) ^^(^T Mi)(*3 ; - Ml) 

^ g (* 2i " M 2 )(*, - Mi) ^ J - Mi) 2 {X lt - ^XX 3i - ^3) 



— g (* „ - M 2 )(*, - Mi) — g (* 3i " M 3 )(* 2i " M 2 ) — ^ (* 3f " M 3 )* 



(4) 



where 'n' stands for a number of elements of a group of color signals. When 'q 1 
pieces of color names are used for color naming, color vector X is to have 'q ! 
pieces of normalized distances. 'S' stands for extensions of a color category's 
distribution in the directions of 'L\ 'a' and 'b\ 

The normalized distance includes the extension of a color category, thus it 
is presumable that a viewer perceives color-vector X to be which color. 

In other words, when the normalized distance between color vector X and 
color category 'red' is smaller than that between color vector X and color 
category 'pink 1 , it might be determined that the viewer perceives color vector X 
to be 'red' rather than 'pink\ 

Euclidean distance in the CIELAB space is used to calculate a distance 
between the gravity of color category T and color-vector X; however, as shown in 
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Fig. 3, color category 'red' has a smaller distribution area than color category 
•pink'. In such a case, because this method does not consider an extension of a 
color category, the color category 'red* would be underestimated if this method is 
directly used for comparing the distances. As a matter of course, a relation 
5 between the distance and the viewers color categorical perception does not 
match. 

On the other hand, the normalized distance absorbs the three- 
dimensional extension difference of color category with the variance covariance 
matrix, the relation between the distance and the viewer's color categorical 

*o 

=J3 10 perception matches well. This is an advantage of using the normalized 

-■=_j 

□ distance. 

s p Next, gravity mapping controller 704 shown in Fig. 7 is detailed. 

01 Fig. 8 illustrates an inside structure of gravity mapping controller 704 in 

accordance with the first embodiment of the present invention. 
15 I* 1 Fig. 8, normalized distance input section 801 for the source device 

outputs a normalized distance. 
O Color-naming database input section 802 for the source device outputs (a) 

each colorimetric value of respective colors 702 for the source device, (b) a 

gravity vector, and (c) a variance covariance matrix all of these are stored in 
20 database 702. Color naming database input section 803 for the output device 

outputs (a) each colorimetric value of respective colors for the output device, (b) 

a gravity vector, and (c) a variance covariance matrix, all of them are stored in 

database 707. 

Gravity-mapping-scaling-coefficient-matrix-calculator 804 calculates a 
25 gravity-mapping-scaling-coefficient-matrix by using the gravity vector and 
variance covariance matrix of database input sections 802 and 803. 

Weighing-coefficient-calculator 805 for the gravity-mapping-scaling- 
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coefficient-matrix calculates a weighing-coefficient which weighs a gravity 
mapping scaling coefficient of each color category based on the normalized 
distance of the source device's color. 

Intermediate normalized distance calculator 806 converts the normalized 
distance of the colorimetric value of the CRT display into an intermediate 
normalized distance by using (a) gravity-mapping-scaling-coefficient-matrix, i.e. 
an output from gravity-mapping-scaling-coefficient-matrix-calculator 804, (b) 
the weighing coefficient for the gravity-mapping-scaling-coefficient-matrix, i.e. 
an output from weighing-coefficient-calculator 805. 

Calculator 806 outputs the intermediate normalized distance to dynamic- 
range-mapping controller 705 shown in Fig. 7 via intermediate normalized 
distance output section 807. 

Next, dynamic-range -mapping controller 705 shown in Fig. 7 is detailed. 

Fig. 9 illustrates an inside structure of a dynamic range mapping 
controller 705 in accordance with the first embodiment of the present invention. 

In Fig. 9, intermediate normalized distance input section 901 outputs the 
intermediate normalized distance. 

Color-naming-database-input-section 902 for the source device outputs (a) 
each colorimetric value of respective colors for the source device, (b) a gravity 
vector, and (c) a variance covariance matrix, all of them are stored in database 
702. Color naming database input section 903 for the output device outputs (a) 
each colorimetric value of respective colors for the output device, (b) a gravity 
vector, and (c) a variance covariance matrix, all of them are stored in database 
707. Source device's color input section 904 outputs the colorimetric value of 
the CRT display's color. 

Dynamic range mapping scaling coefficient matrix calculator 905 
calculates a dynamic-range-mapping-scaling-coefficient-matrix by using (a) the 
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gravity vector, (b) the variance covariance matrix, both supplied from color- 
naming-database input section 902 for the source device, and (c) the gravity- 
vector, (d) the variance covariance matrix, both supplied from color-naming- 
database input section 903 for the output device. 
5 On the other hand, dynamic range mapping restraint coefficient 

calculator 906 calculates a dynamic range mapping restraint coefficient by using 
(a) the gravity vector, (b) the variance covariance matrix of input sections 902 
and 903, and (c) the colorimetric value of the CRT display's color. 

Weighing coefficient calculator 907 for the dynamic-range -mapping- 
*0 10 scaling-coefficient-matrix calculates a weighing coefficient for the dynamic- 
p range-mapping-scaling-coefficient-matrix by using the colorimetric value of the 

CRT display's color. 

Target normalized distance calculator 908 receives the intermediate 
normalized distance from input section 901, and converts the received distance 
1™ 15 to a target normalized distance by using the dynamic-range-mapping-scaling- 
Im coefficient-matrix weighing coefficient. 

O Target normalized distance output section 909 outputs the target 

normalized distance to target-normalized-distance-searching-section 706 shown 
in Fig. 7. 

20 Away of producing database 702 for the source device or database 707 for 

the output device is described with reference to Fig. 10. 

Fig. 10 shows a color-naming experimental system. 

In Fig. 10, on the CRT display or on a paper 1001 of the printer, at least 
one color patch 1002 for test stimulus is rendered, the colorimetric values of 
25 color patch 1002 have been measured. Viewer 1003 watches this color patch 
and names the colors with a plurality of color names. To be more specific, the 
color names to be selected by the viewer are rendered on the CRT display or on 
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the paper 1001 by the printer. The selected color names of color patch 1002 are 
input with mouse 1004 into calculator 1005. Instead of the mouse, voice 
recognition can be used for the input. A number of viewers is not necessarily 
one, but it may be plural. 
5 When color patch 1002 is displayed on CRT display 1001, calculator 1005 

stores the colorimetric values of color patch 1002 — corresponding to the color 
names input with mouse 1004 by the viewer — into color-naming-database 702 
for the source device. Instead of the colorimetric values, device-drive-levels can 
be stored. 

10 When color patch 1002 is printed on a paper, calculator 1005 stores the 

colorimetric values of color patch 1002 — corresponding to the color names input 
with mouse 1004 by the viewer — into color-naming-database 707 for the output 
device. Instead of the colorimetric values, device-drive-levels can be stored. 

It is effective to select a color name for the color naming responsive to an 

15 objective of color reproduction. For instance, when color categorical 
characteristics should be described regarding an entire color space of an image 
device, the basic 11 color category names found by Berlin and Kay is 
recommended. (Re: basic color category names, reference document: Basic Color 
Terms, "Their Universality and Evolution", Univ. of California Press, Berkley, 

20 1969) 

Nearly 100 languages were studied, and it has been found that any 
developed languages commonly share these 11 basic categorical color names. 

Mr. Crawford summarizes a definition of the basic color names as follows: 
(reference document: "Defining basic color terms", Anthropol Linguist, p24, 
25 p338-343, 1982) 

1. The basic color names must be included in every language. 

2. The basic color names are independent of users as well as cases, 
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and they must be used in a stable manner. 

3. The meaning of the basic color name is not included in other words. 

4. The basic color names must be used for a specific object. 

Mr. Uchikawa introduces 11 basic color categorical names in Japanese. 
5 (reference document: "Perception and Memorizing of categorical colors", the 7th 
color engineering conference, 1990) If such basic categorical color names are 
used, color names can be selected with respect to the entire color space of an 
image device without overs and shorts. 

One of critical factors for determining the color reproducibility is 
10 'memorized colors'. When mapping control of the memorized colors is used as 
an important tool, it is effective to use the memorized colors as color names. 
Setting color names indicates geometrically that a gravity is set in a color space. 
Therefore, when the memorized color such as 'sky blue' or 'grass green' is used 
as a color name, it is expected that the color reproducibility is improved because 
15 mapping points are controlled based on the memorized color, which is critical for 
color reproduction. 

When the viewer uses a color name— perceived at watching a test image 
— for the color naming, an appropriate color-information can be exchanged free 
from loss comparing with a case where an entire gamut of the device is collected, 
20 because the mapping can be controlled to place some weight on an image gamut. 

Categorical basic colors, memorized colors and the colors perceived by a 
viewer are respectively used for the color naming. Effectiveness of the 
respective cases is discussed above. Proper color categories other than the 
above also can be set. 

25 The mapping process to the printer's color space with preserving the color 

categorical characteristics of the CRT display's color space is detailed 
hereinafter with reference to Fig. 11. 
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Fig. 11 shows a procedure of mapping into printer's space with preserving 
color categorical characteristics of the color space of a CRT display. 

In this procedure, a number of color categories is 9 for the convenience. 
(White, black and gray are deleted because they do not have chroma.) 
5 In Fig. 11, the procedure starts by inputting a drive signal from CRT 

display driving unit 103 to the CRT display (source device) in step 1101, then the 
drive signal is converted to the colorimetric value of an emitting color in step 
1102. 

Input CRT display's gravity vector U &i , and variance covariance matrix 
jS 10 E si , in step 1103, both are found by color naming, then further input printer's 
O gravity vector fi s ii and variance covariance matrix 2 si , in step 1104. 

5 p Calculate normalized distances with respect to the gravity vectors of the 

m nine c °l° r categories with equation (1) using the CRT display's colorimetric 

u value, gravity vector £t Sii , and variance covariance matrix S si in step 1105. 

7 15 The nine normalized distances found in step 1105 are converted to 

intermediate normalized distances based on a positional relation between the 
O CRT display's vector and the printer's vector in step 1106. Then the 

intermediate normalized distances are converted to target normalized distance 
in step 1107 so that optimum points can be mapped by a dynamic range 
20 mapping operator based on mapping information on a surface control point 
provided on a gamut surface or near to the gamut surface. 

Through the process in step 1107, the colors on the CRT display (source 
device) are mapped on the printer (output device) as nine normalized distances. 
Then, search the color space of the printer for the colorimetric values 
25 having the target normalized values in step 1108. 

Printer driving unit 105 calculates a drive level using a printer profile in 
step 1109 in order to print the colorimetric values searched, then outputs the 
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drive signal to the printer in step 1110. 

Next, the CRT display's normalized distance vector D s is converted to 
target normalized distance vector D t . This process is described with reference 
to equations. Equation (5) converts vector D s to vector D t with gravity mapping 
operator 'V and dynamic mapping operator 'R\ 
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With gravity mapping operator 'V, vector D s is converted to distance 
vector D btt , namely the intermediate normalized vector. 

The elements V 1( V 2 , ,V 9 of the gravity mapping operator 'V is 

extracted from equation (6). 
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where $ is referred to as a gravity mapping scaling coefficient matrix, and this 

is extracted from equation (7). Also [w vl , w v2) , w^]' are the weighing 

vectors for the gravity mapping scaling coefficients, and each element is 
extracted from equation (8). 



<t> = 



V 21 V 22 



V 91 V 92 



19 



29 



99 



jm 



S.jm 

0 



(7) 



(l*m) 
(Z=m) 



/,m=l,2,--,9 



1 

W ■ - DsJ 

v > ; 9 1 

if D s> , = 0 then w v> , = 1, w v> k m t = 0 (8) 
*/ &s,k* i =0 then w Vik „ = l,w v> , „ f 

Gravity mapping scaling coefficient matrix $ is an operator for mapping 
the CRT display's gravity to the gravity of the same color category of the printer. 
Geometrically, this matrix rearranges the positional relations between nine 
gravities of the CRT display to those of the printer. 

Element v lm of gravity mapping scaling coefficient matrix <1> is a scaling 
coefficient converting normalized distance D s lm to another normalized distance 
i m is the normalized distance viewed from the gravity of color category 
T toward the gravity of color category 'm* in the color space of the CRT display, 
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and D., lm is another normalized distance viewed from the gravity of color 
category T toward the gravity of color category 'm' in the color space of the 
printer. Gravity mapping scaling coefficient matrix <I> includes all scaling 
coefficients corresponding to every possible combination of the nine color 
5 categories, therefore, all the gravity vectors in the CRT display's color space are 
mapped at the points of gravity of the same color category in the printer's color 

space. However, as equation (7) shows, the scaling coefficient (v u , v mm ) for 

itself is zero (0). 

Gravity mapping for other than the gravity is carried out by interporating 
10 the mapping information, the gravity has, responsive to the positional relation 
between color vector X and gravity vector PL . 

In equation (6), column vector elements [v H , v 2i , , v 9 J of gravity 

mapping scaling coefficient matrix $ are the scaling coefficients realizing the 
mapping between gravities in color category ? i', and they are weighed by 
15 weighing coefficient w vi . As shown in equation (8), weighing coefficient w vi is 
defined by B/A, where 'A' represents the sum of inverses of the normalized 
distances between the color vector X and the gravities of nine color categories, 
and 'B 1 represents the inverse of the normalized distances between the color 
vector X and the gravities of color category T. Accordingly, the gravity of the 
20 color category nearer to color vector X has more important scaling coefficient, 
and the gravity of the color category farther from color vector X has less 
important scaling coefficient. 

As discussed above, gravity mapping operator 'V 1 is extracted from 
equations (6), (7) and (8), and operator 'V converts the CRT display's normalized 
25 distance D s to intermediate normalized distance D b t . 

In equation (5), dynamic range mapping operator 'R' converts 
intermediate normalized distance vector D b t into target normalized distance 
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vector D t . The function of operator 'R* is to incorporate the difference in gamut 
shapes between the CRT -display (source device) and the printer (output device) 
into the intermediate normalized distance vector D bt which only underwent the 
mapping between the gravities by operator 'V 1 . This process corresponds to a 
5 dynamic range mapping. 

Operator 'R', as shown in equation (5), comprises dynamic range 
correction coefficient U i ~ making the dynamic range of the source device 
correspond to that of the output device — and restraint coefficient F which 

^ restrains the work of coefficient Uj. 

jfjj 10 First, coefficient Uj is detailed hereinafter. 

™ In. order to design dynamic range correction coefficient JJ i9 a surface 

control point must be defined. The surface control point is disposed on a 

m surface of device's gamut or near to the surface of the gamut. The point 

functions as a mapping target of the dynamic range mapping. 

j*i 15 I n this embodiment, for an easy description, eight surface control points 

jSj are defined; however, a number of the control points can be set arbitrarily. 

™ The eight surface control points are defined in digital counts of the device 

as follows: the primary color: red R:(d r , d g , d b ) = (255, 0, 0), green G: (d r , d g , d b ) = 
(0, 255, 0), and blue B: (d r , d g , d b ) = (0, 0, 255). Secondary colors: yellow Y: (d r , d g , 
20 d b ) = (255, 255, 0), magenta M: (d r , d g , d b ) = (255, 0, 255), cyan C: (d r , d g , d b ) = (0, 
255, 255), white W: (d r , d g , d b ) = (255, 255, 255), black B: (d r , d g , d b ) = (0, 0, 0). 
These eight points are positioned at respective vertexes of the device's eight 
gamuts, and on the surface of the gamuts. These surface control points are set 
both on the CRT display and the printer, and the mapping is carried out 
25 between the control points of the same color. This is done by dynamic range 
correction coefficients U { . Therefore, coefficient 1^ functions as mapping again 
the surface control points of the CRT display onto the surface control points of 
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the printer. The control points of the CRT display have been once mapped on 
the CRT display by gravity mapping operator *V\ 
Coefficient XJ { is extracted from equation (9). 
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where "Sf is referred to as a dynamic range mapping scaling coefficient matrix 

and extracted from equation (10). Weighing vector [w r l , w r2 , , w r9 ] for 

the dynamic range mapping scaling coefficient is extracted from equation (11). 
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where ! p' represents a kind of the color categories and 'q' represents a kind of the 
surface control points. 

Dynamic range mapping scaling coefficient matrix NP" is an operator that 
maps the surface control points of the CRT display onto the surface control 
points of the same color of the printer. This operator works geometrically 
equivalent to gravity mapping operator V that performs gravity-to- gravity 
mapping. In other words, matrix Y rearranges the positional relations 
between eight surface control points of the CRT display to those of the printer. 
Element i/> pq of matrix"** is a scaling coefficient that converts CRT display's 
intermediate normalized distance D dabpq into printer's intermediate normalized 
distance D d a c pq . Distance D d a b pq extends between the surface control point of 
color ! q' and the gravity of color category 'p' in the color space of the CRT display, 
and distance D dacpq extends between the surface control point of color 'q' and the 
gravity of color category f p' in the color space of the printer. Matrix includes 
all the scaling coefficients with respect to every combination of the nine color 
categories and the eight surface control points. However as equation (11) 
shows the scaling coefficient (p = q) for itself is zero (0). 

Dynamic range mapping for other than the surface control points is 
carried out by interporating the mapping information, the surface control point 
has, responsive to the positional relation between color vector X and surface 

control point. In equation (9), column vector elements u,^^ >^sJ of 

dynamic range mapping scaling coefficient matrix are the coefficients which 
scales the intermediate normalized distance from gravity points in color 
category T to the surface control points, and they are weighed by weighing 
coefficient w r i . As shown in equation (11), in the color space of the CRT display 
weighing coefficient w ri is defined by D/C, where f C represents the sum of 
inverses of the Euclidean distances in CIELAB space between the color vector X 
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and the eight surface control points of color vector X, and 'D* represents the 
inverse of the Euclidean distances between the color vector X and the surface 
points of color category ! i\ Accordingly, the surface control point of the color 
category nearer to color vector X has more important scaling coefficient, and the 
5 surface control point of the color category farther from color vector X has less 
important scaling coefficient. 

Then another element, i.e. restraint coefficient T\ of dynamic range 
mapping operator 'R' is detailed hereinafter. 

Dynamic range correction coefficient U { maps the surface control points of 

10 the CRT display (source device) onto the surface control points of the printer 
(output device). However, in this case, CRT display's gravities that were 
mapped onto the gravities of the printer by gravity mapping operator 'V are 
again mapped by operator 'R', thereby moving to other places than the printer's 
gravities. If the design policy of mapping the gravities between the CRT 

15 display (source device) and the printer (output device) as well as mapping the 
surface control points therebetween is carried out, the color vector once mapped 
to the printer's gravity by operator 'V must be controlled not to be affected by 
the mapping by operator 'R ? , so that the color vector does not move from the 
gravity. 

20 If equation (5) is reconsidered, intermediate normalized distance vector 

D b t converted by operator 'V directly becomes target normalized distance vector 
D t , when color vector X agrees with the gravity of the CRT display (source 
device). 

Accordingly, correction coefficient JJ i should not work at all when color 
25 vector X agrees with the gravity of the CRT display (source device), but it should 
work to the maximum extent when color vector X agrees with the surface control 
point. In order to make coefficient \J i work in such a way, the effect of 
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coefficient U i must be weighed considering the positional relation between color 
vector X and the gravity point of the CRT display (source device). Restraint 
coefficient 'F' will do this job, and coefficient 'F is extracted from equation (12). 
F=Wm.ifi+ w mt2 f 2 + • • • + w mt9 f 9 (12) 

where fjis given by equation (13) and w mi is given by equation (14). 
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20 where E s m i represents Euclidean distance in the CIELAB space between gravity 
vector fi i of color category 'iV and color vector X. E s ai represents Euclidean 
distance in the CIELAB space between surface control point *i' and color vector X. 
Both of E sm i and E sai indicate the distance in the source color space (color space 
of the CRT display). According to equation (13), %' works as binding force to 

25 correction-coefficient Ui which is given by the gravity vector of color category V 
to the eight surface control points. When E s m i = 0, namely color vector X agrees 
with gravity vector u i3 then £ becomes 1 and the maximum binding force is 
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expected. On the contrary, when E s a i = 0, namely, the test color agrees with the 
surface control point, then f { becomes zero (0) and no binding force is expected. 

On the other hand, in equation (14), D sm i represents the normalized 
distance between gravity vector tt { of color category V and color vector X. This 
distance is limited in the source color space. When color vector X agrees with 
gravity vector u it w m i becomes 1 (one) and weighing coefficient w mk!=i becomes 
zero (0). Therefore, w m (i functions as a weighing coefficient to binding force f { 
determined according to the position of color vector X. Respective binding 
forces %' of the nine gravity vectors are weighed by w m ti and the results are 
reflected to restraint coefficient ! F\ 

Through the process discussed above, the input colorimetric value of the 
CRT display is converted to the target normalized distance, and the mapping 
based on the normalized distance is completed. 

Step 1108 shown in Fig. 11, i.e. the process to find a colorimetric value 
having the target normalized distance out of the values in the color space of the 
printer is detailed hereinafter. 

Fig. 12 shows a procedure for searching the color space of the printer 
(output device) for a colorimetric value having the gravity vector and the target 
normalized distance. 

Step 1108 begins with inputting the target normalized distance (step 
12001). 

Then input the gravity vector and a variance covariance matrix of the 
printer (output device) in step 12002, and set a candidate color for outputting in 
step 12003. 

Calculate the normalized distance between the candidate color and the 
gravity vector of the printer in step 12004. 

Give the normalized distance to a function of the candidate color for 
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outputting, i.e. equation (15), in step 12005. Then compare the function with a 
predetermined threshold value in step 12006. 

When the function shows a greater value than the threshold value, renew 
the candidate color in step 12008, and return to step 12004. On the other hand, 
the function shows a smaller value than the threshold value, end the search, 
and send the candidate color to the output device as an output colorimetric value, 
then complete the process in step 12007. 

f(p) i( D -'- D ")' , { °->-D'A ' , , (D.>-D,A 2 

where D pi represents has the normalized distance between the candidate color 
and the gravities of color category T. 

Meanwhile, the CMS needs to change its form following an objective of 
color reproduction. For instance, when colorimetric matching is the objective of 
the color reproduction, the CMS is optimized based on the reference that a 
source color agrees with an output color in terms of CIEXYZ value or CIELAB 
value. When a matching of color appearance is the objective of the color 
reproduction, the reference should be the agreement of the attributions of the 
source color and output color, where the attribution includes brightness and 
colorfulness in absolute term, and hue, lightness, chroma in relative term. 

In a desirable color production, an agreement of memorized colors with 
reproduced color should be a key, and the color should be reproduced according 
to viewer's perceptive color. 

The first embodiment of the present invention can be selectively applied 
to a variety of color reproductions by using gravity-control points formed by 
gravity vectors as well as the surface control points disposed on gamut surfaces 
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or near to the gamut surfaces. 

In other words, when a vector agreeing with the source device is 
colorimetrically used at a control point of the output device, colorimetric 
matching is preserved between the control points of the source and output 
5 devices, so that the mapping intermingled with colorimetrical matching and 
preserving of categorical characteristics for gamut matching is realized. 

In a section near to gray axis in the gamut, colorimetrical matching is 
realized, and on the gamut surface, a non-linear process due to gamut 
compression is realized. Colors have been excellently reproduced in that way 
10 in general. However, in the present invention, gravity control points in the 
gamut form colorimetrical-mapping-pair, and the surface control points form 
color-perception-mapping-pair (e.g. color appearance mapping), then an ideal 
mapping process is realized. 

If differences in color appearances are minimized on every control-points 
15 pair, a gamut mapping — considering chromatic adaptation due to difference in 
lighting — can be designed, and a gamut mapping method to be optimized color 
categorically to arbitrary viewing environment can be designed. 

When the first embodiment does not accompany the gamut compression 
in mapping, it works as a color-information-exchanging-device that can realize a 
20 color appearance matching, i.e. predicting a corresponding color. 

Specifying the color — of which color and appearance agreed with each 
other on a source device — on an output device requires a process where a viewer 
selects the color out of plural given colors. 

When the color cannot be specified on the output device due to different 
25 gamut of the source and output devices, a process is required where a color 
showing a minimum difference of appearance is searched. 

Fig. 13(a) shows a cusp of gamut of a source device and that of an output 
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device. The horizontal axis represents hue, and vertical axis does lightness. 

If the colors between cusp 1301 of the source device's gamut and cusp 
1302 of the output device's gamut are used as test colors, the process discussed 
above can be finished in a short time. 

Fig. 13(b) shows gamut boundaries 1303 and 1304 of the source device 
and output device respectively at hue angle of 200 degree, where the boundaries 
are indicated by a plane of chroma C*~ lightness L*. 

When a surface control point of the source device is disposed at cusp 1305 
where the maximum chroma is rendered on gamut boundary 1303, the 
corresponding color having a minimum difference in appearance from the 
control point seems to be around cusp 1306 where the maximum chroma is 
rendered on gamut boundary 1304 of the output device. 

If cusp 1306 per se is not selected, it is natural to think that a section 
nearer to the source device's cusp 1305, i.e. brighter section, can be selected 
rather than a section darker than cusp 1306. Accordingly, the color where the 
difference in color appearance is minimized is possibly between cusp 1301 and 
cusp 1302. Therefore, if the test colors are selected out of this region, the color 
which shows minimum difference in color appearance or agreement in color 
appearance can be searched within a short time. 

As discussed above, the first embodiment proves that the gamut mapping, 
where color attributes and color reproducibility are quantatively understood, 
can be realized, and also a practical corresponding-color-prediction, absorbing 
the difference in viewing conditions, can be realized. The procedure discussed 
in the first embodiment is this: at least one viewer names colors, then the colors 
are grouped based on colorimetric values or device driving level of the test colors 
into color categories. Each color category is normalized by using statistical 
characteristics such as variance, so that normalized color space is provided. In 
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this color space, control points of the same color name are mapped, so that a 
source device and an output device exchange color information. Through these 
procedures, the two features discussed above can be realized. 

The first embodiment also proves that it can minimize the difference in 
5 respective impressions of an entire image on the source device and the output 
device through the following procedure: Classify the CIELAB color space of the 
source and output devices into a plurality of color categories, and p r °cl u ce 
normalized color spaces for the respective devices by using (a) a gravity vector 
q calculated in each color category and (b) variance covariance matrixes 

10 expressing extension from the gravity points. Then map the control points of 
^ the same color in the normalized color spaces of the respective devices, thereby 

;=F exchanging the color information by preserving the agreement of the color 

\7i categories in the two normalized color spaces as well as the sizes of color 

perception quantity such as chroma and the like in the two normalized color 

H ' 
1^ 15 spaces. 

IjJ 

Exemplary Embodiment 2 

The color information exchanging device in accordance with the second 

exemplary embodiment is demonstrated hereinafter with reference to Fig. 1 and 
20 Fig. 14. This is the case where a CRT display and a print by a printer differ 

greatly in a reference white color, i.e., mapping an color appearance is performed 

in this embodiment. 

Fig. 14 illustrates an inner structure of the color information exchanging 

device used in the second embodiment. 
25 In Fig. 1, CRT display driving unit 103 outputs the colorimetric values of 

the colors displayed on CRT display 102. These colorimetric values are 

converted to lightness, chroma, and hue, which are three attributes of human 
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visual sight to the colors, by polar coordinate converter 2104 shown in Fig. 14, 
then the converted values are supplied to two systems. 

One of the systems is gamut mapping per se. In this system, gamut 
mapping undergoes pre-mapping section 2109 and main-mapping section 2110. 

The other system is used for finding a gamut-shape-control-coefficient. 
In this system, the converted values undergo max.-chroma-color-setting-section 
2105 where the max. chroma of the colors on the CRT display is set on gamut 
boundary, pre-mapping section 2106 where the colors on the gamut boundary 
are mapped on gamut boundary of a print by the printer, main-mapping section 

2107 for the max. chroma color, and gamut-shape-control-coefficient-calculator 

2108 where an adjusting coefficient — allowing the gamut boundary of the CRT 
display to be mapped on the gamut boundary of the print — namely, the gamut 
shape control coefficient is calculated. This resultant coefficient is supplied to 
main-mapping section 2110 and used for main mapping. An output from 
main-mapping section 2110 is converted from lightness, chroma and hue in the 
polar coordinate into colorimetric values in an orthogonal-coordinate-converter 
2111, then the resultant values are supplied to printer driving unit 105 shown in 
Fig. 1. 

The processes in the pre-mapping and main-mapping are described in Fig. 
15. Fig. 15 illustrates processes of the main-mapping and pre-mapping when 
lightness mapping from a CRT display to a printer is performed in accordance 
with the second exemplary embodiment of the present invention. Inside of 
gamut boundary 2201 is a gamut of the CRT display. Insides of gamut 
boundaries 2202 and 2205 are printer's gamut at hue H s . Black point 2203 
represents a color of the CRT display, black points 2204 and 2207 are pre- 
mapped colors at hue H s , inside of gamut boundary 2206 is the printer's gamut 
at hue H d , and black point 2208 represents a color by the printer at hue H d . 
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At pre-mapping section 2109, gamut 2201 of the CRT display is once pre- 
mapped into printers gamut 2202, then p re-mapped color 2207 at hue H s is 
mapped again at the position of black point 2208 in hue H d so that the color 
category of color 2203 of the CRT display can agree with the color category of 
5 printer color 2208. 

Pre-mapping section 2109 has a function of mapping along a geometrical 
construction of a gamut, and main-mapping section 2110 has a function of 
maximizing the color reproducibility. However, when the main-mapping is 
performed, the geometrical construction of the printer's gamut must be obtained, 

10 so that the colors once delivered into the printer's gamut by the pre-mapping 
will not go out from the printer's gamut. Gamut-shape-control-coefficient 
calculator 2108 has a function of supplying information of gamut shape to 
main-mapping section 2110, and allows the printer colors determined by the 
main mapping to be situated in the printer's gamut. 

15 The concept of the pre-mapping is described hereinafter with reference to 

Fig. 16. Fig. 16 is a cross section of color space defined by lightness L*, chroma 
C and hue H taken on line L*~C* at hue H s of an emitting color of the CRT 
display. The defined color space has been converted from the CIELAB by polar 
coordinate conversion. The left side on Fig. 16 is a cross section of CRT 

20 display's gamut, and inside of gamut boundary 2301 is a gamut of the CRT 
display. The right side is a cross section of printer's gamut. Inside of gamut 
boundary 2302 is printer's gamut. 

Pre-mapping includes the following steps: first, determine white point 
2303 on the lightness axis and white point 2304 on the gamut boundary, then 

25 put both white points 2303 and 2304 at a chromatic saturated point of the color 
on the CRT display, and finally, determine a pre-mapping point (represented by 
a black point). This is shown by equation (16). 
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The lightness L* u g of white point 2303 and the lightness L* u surface of white 
point 2304 are found in this way. 



L*' 




C 




H 


H 



(l-C,)L\,s+L\, 

CrrC^d, max 

Hs 



surface 



(16) 



where C* u , dmax represents the maximum chroma at hue H s , light L UtSurface . C* r is 
given by equation (17). 
10 C = C* s /C* s , max (17) 

where C* s represents chroma of a color of the CRT display, C* s max represents the 
chroma of the color having the maximum chroma in the gamut of the CRT 
display, the gamut having the same luminance and the same hue as the color by 
15 the CRT display. 

Lightness L* ug is given by equation (18), i.e. lightness L* s of the CRT 
display's color 2305 is linearly compressed. 



20 



s 



L* 



d, min 



[L dt max L 



d t min 



in/ T * 



(18) 



T* -T * 

J^f 5, max ^ s, 
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where L* s max represents the maximum luminance of the CRT display's gamut, 
L* s , m in represents the minimum luminance thereof, L* d max does the maximum 
luminance of printer's gamut and L* d min does the minimum luminance thereof. 

On the other hand, Lightness L* u surface is extracted by supplying lightness 
L* s of color 2305 to equation (19). 
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7*7* 

T * - t * +fr* _r * ^ ^ 5 5, mm / r * r * \ 

u F jur/acff 1^ rf, min U-f d, cusp , mapped ^ rf, min/~"T^ ~ s ^ S ,CUSp , mapped } ^ 

-Li s, max ^ J, min 
7 * _7 * 

T * +(T* 7* \ * * min (7* 7* \ 

^ d, min rf, max ^ < min / + # I J~J s JU S,CUSp, mapped 

<Li Lj s, min ' 



where L\ ^ represents the luminance of the CRT display's cusp, and L* s cusp , 
mapped represents the lightness of CRT display's cusp mapped on the gamut 
surface of the printer. The cusp is a sharp point of the gamut having the 

10 maximum chroma in a hue, and point 2306 in Fig. 16 is a CRT display's cusp at 
hue H s , and point 2307 is a printer's cusp at hue H s . 

Lightness L* s cusp mapped of the CRT display's cusp mapped on the gamut 
surface of the printer is calculated in this way: Fig. 17 illustrates cusps of the 
CRT display and the printer. Fig. 17(a) is a projection drawing of respective 

15 cusps of a CRT display and a printer of CIELAB space projected to a plane of a* 
— b*. Line 2401 represents the cusp of the CRT display and line 2402 does that 
of the printer. Fig. 17(b) is a projection drawing of respective cusps of a CRT 
display and a printer of CIELAB space converted to a polar coordinate projected 
to a plane of L* — H (lightness — hue). Line 2403 represents the cusp of the CRT 

20 display and line 2404 does that of the printer. As shown in Fig. 16, since the 
cusp has a high chroma in a middle luminance, it draws viewer's attention. 

From the viewpoint of color reproducibility, controlling the cusp is a key 
for designing a gamut mapping. 

Fig. 18 is a schematic diagram of Fig. 17. 

25 In Fig. 18(a), line 2501 represents the CRT display's cusp, and line 2502 

does that of the printer. Line 2503 represents the CRT display's cusp mapped 
on the gamut surface of the printer. Black point 2504 is a color displayed on 
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the CRT display, white point 2505 represents a mapping color of point 2504 
mapped on the surface of the printer's gamut by equation (19). Fig. 18(a) 
shows a case where the lightness L s of the color by the CRT display is lower than 
the lightness L* 8t ^ mapped of the CRT display's cusp mapped on the surface of the 
5 printer's gamut. On the other hand, in Fig. 18(b), lines 2506, 2507 and 2508 
represent the CRT display's cusp, printer's cusp and the CRT display's cusp 
mapped on the surface of printer's gamut, respectively. Black point 2509 
represents a color displayed by the CRT display, white point 2510 represents a 
mapping color of black point 2509 mapped on the surface of the printer's gamut 
10 by equation (19). Fig. 18(b) shows a case where the lightness L* s of the color by 
the CRT display is higher than the lightness L* St CUSPi mapped of the CRT display's 
cusp mapped on the surface of the printer's gamut. The cusp data (2501, 2502, 
2503) shown in Fig. 18(a) is identical to the cusp data (2506, 2507, 2508) shown 
in Fig. 18(b). 

15 When the lightness L* s of the color by the CRT display is lower than the 

lightness L* s CUSP) mapped of the CRT display's cusp mapped on the surface of the 
printer's gamut, i.e. in the case of Fig. 18(a), 'C is extracted from equation (20) 
and L* UjSurface is extracted from equation (21). 

(A/B)=(C/D) (20) 

20 

L* usurface = L* dmin +C (21) 

This case agrees with a case of L*s □ L* s CUSP( mapped in equation (19). On the 
other hand, when the lightness L* s of the color by the CRT display is higher than 
25 the lightness L* s> cusPt mapped of the CRT display's cusp mapped on the surface of 
the printer's gamut, i.e. in the case of Fig. 18(b), 'G' is extracted from equation 
(22), and L* St surface is extracted from equation (23). This case agrees with a case 
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ofIA>L* 



s, cusp, mapped 



in equation (19). 



(E/F)=(G/H) 



(22) 
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u, surface' 
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+G 



(23) 



Next, finding L* s cusp mapped (the lightness of the CRT display's cusp mapped 
on the printer's gamut) in equation (19) is described hereinafter. 

As Fig. 17 shows, cusp 2403 of the CRT display has a lightness different 
from cusp 2404 of the printer. This phenomenon can be seen between many 
image devices. 

The present invention thus considers a color category distribution on a 
gamut surface, then maps the cusp of the CRT display onto the gamut surface of 
the printer. 

Fig. 19(a) shows a color category distribution of color names that are 
named by a viewer with respect to the colors on the gamut surface of the CRT 
display. 

Fig. 19(b) shows a color category distribution of color names that are 
named by a viewer with respect to the colors on the gamut surface of the printer. 
In Fig. 19, the relation between symbols and colors is this: 



white triangle A : red 



black plus 



+ : brown 



black cross X : pink 



black ball 0 : orange 

white wedge O : yellow 
black triangle A : green 
white square □ : blue 



black wedge 



♦ : violet 
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In Fig. 19, the common feature of the CRT display (Fig. 19(a)) and the 
printer (Fig. 19(b)) is that each color has a different spread of distribution in the 
lightness direction. For instance, green (A) and blue (□) are spread widely in 
the lightness direction; however, yellow(<C>) and red (A) are spread narrowly. 
Therefore, when the colors on CRT display's gamut surface are mapped onto 
printer's gamut surface, narrow spread colors such as yellow (O) and red (A) 
must be mapped within a limited rang so that the color category on the CRT 
display's gamut surface can agree with the counterpart on the printer's gamut 
surface. For instance, the yellow (O) spread in the range of 60 < L < 95 on the 
CRT display should be mapped in the range of 80 < L < 95 so that the color 
category appears yellow on a print by the printer. When the yellow color of L = 
70 on the CRT display is mapped in the printer's space with preserving the 
lightness, the color category turns to brown (+). If the mapping is viewed from 
the standpoint of chroma, the CRT display's cusp must be mapped into the 
printer's cusp in order to minimize the lowering of chroma. It can be concluded 
that the two objectives, i.e. agreement between color categories and preserving 
of chroma, are balanced appropriately in determining mapping points. 

Thus the present invention allows the mapping in the following way: 
Among the cusps of a source device, the colors widely spread in the lightness 
direction are mapped near to the cusps of the source device, and the colors 
narrowly spread are mapped close to the cusps of an output device. For 
instance, yellow (O) and red (A) are spread narrowly in the lightness direction, 
thus among the CRT display's cusps, cusps belonging to yellow (O) and red (A) 
are mapped on the printer's cusps. On the other hand, since green (A) and 
blue (□) are spread widely in the lightness direction, among the CRT display's 
cusps, the cusps belonging to green (A) and blue (□) are mapped on the cusps of 
the CRT display. 
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In order to realize this idea, cusp-lightness-difference-correction- 
coefficient !f is introduced. The lightness L s cusp mapped of the CRT display's cusps 
mapped on the printer's gamut is extracted from equation (24); 

s.cusp.mapped ~~ L d,cusp *^~te (24) 

5 e = L\ cusp - L\ cusp s 0 

where L* d ^ represents luminance of the printer's cusp. 

The lightness L* s ^ mapped of the CRT display's cusps mapped on the 
printer's gamut is closer to CRT display's cusp at greater 't' (cusp-lightness- 
10 difference -correction-coefficient), on the other hand it is closer to the printer's 
cusp at smaller ! t\ 

According to the present invention, for the colors spread widely in the 
lightness direction, cusp-lightness-difference-correction-coefficient 't' increases, 
and for the colors spread narrowly, coefficient 't' decreases. Cusp-lightness- 
15 difference-correction-coefficient 't' is extracted from equation (25). 

t =P „ H >^- H > + Ps , l + 1 H °- H > (i = 1,2,...,6) (25) 
H l + l — H x H i + t —Hi 

where H t represents a hue angle of the 'i'th divided hue, and P St i, P 8 i+ , are 
20 extracted from equation (26). 
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Ps,j = S i LL - f (j =1,2,...,7, Ps,7 = Ps,l) 



LL, max 



^•^H^fc)" 1 ' (26) 



n -1 (A: =1,2,. ..,6) 
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where S Lihk represents a standard deviation of lightness component L k of color 
category 'k\ max( ) represents the maximum value of an element specified in 
parenthesis ( ), and min( ) represents the minimum value of an element 
specified in ( ), where p s 7 = p s { . 

Standard deviation S LL k of lightness component of color category ! k\ is 
supplied from the database 905 of lightness component dispersion. This data- 
base is described later. 

Lightness component variance database 905 is produced in the following 
way: The variances of each color category are calculated with a plurality of 
colorimetric values grouped by viewers through color naming shown in Fig. 10. 
Variance database 905 can be produced in the same way as color naming 
database 707 for the output device by using the system shown in Fig. 10. The 
detailed description is thus omitted here. 

In Fig. 10, viewer 1003 watches color patch 1002 having at least one test 
stimulus of which colorimetric value has been measured on the CRT display or 
the printer. The viewer names the color with a plurality of color names. To be 
more specific, color names of color patch 1002 — answered by the viewer and 
rendered on the CRT display or a print by the printer — are fed into calculator 
1005 with mouse 1004. When color patch is displayed on the CRT display, 
calculator 1005 stores the colorimetric values of color patch 1002 or a device- 
drive-level into lightness-component-variance database 905 responsive to the 
color names input through mouse 1004. 

Variance covariance matrix T, { shown in equation (27) is calculated using 
the colorimetric values of each color stored in database 905. 
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For the color names, categorical basic color names shown in table 1 are used. 

TABLE 1 





_.T T A O 


sL * sa 


sL * sb 


Category 


sL * sa 


saa A 2 


sa * s 




sL * sb 


sL * sb 


sbb A 2 


1 


547.14 


-5.32 


72.39 


Achromatic 


-5.32 


42.50 


-29.22 




72.39 


-29.22 


115.28 


2 


84.91 


92.32 


66.39 


Red 


92.32 


150.97 


83.74 




66.39 


83.74 


171.29 


3 


289.98 


-54.61 


125.78 


Brown 


-54.61 


120.48 


-105.61 




125.78 


-105.61 


346.42 


4 


173.26 


-98.38 


7.37 


Pink 


-98.38 


253.98 


-37.42 




7.37 


-37.42 


83.36 


5 


197.10 


-62.18 


32.23 


Orange 


-62.18 


226.65 


-54.15 




32.23 


-54.15 


320.82 


6 


13.04 


-3.02 


-11.40 


Yellow 


-3.02 


23.64 


-50.15 




-11.40 


-50.15 


340.27 


7 


202.48 


52.16 


18.75 


Green 


52.16 


219.84 


22.72 




18.75 


22.72 


372.51 


8 


335.67 


-191.86 


126.79 


Blie 


-191.86 


295.46 


-129.90 




126.79 


-129.90 


164.13 


9 


223.89 


-97.20 


71.21 


Purple 


-97.20 


223.21 


-45.20 




71.21 


-45.20 


102.39 
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Table 1 shows an example of variance covariance matrix of the printer. In 
order to make every colorimetric value be a positive value, [L*, a*, b*] 1 is 
converted to [a , j3 , 7 ] e with equation (28), and variance covariance matrix Z t 
is found by using [a, j3 , 7 ]\ 



"a" 




L *" 




"78" 






a* 


+ 


128 


y 




b* 
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Lightness dispersion is represented by S LL 2 at the first row and the first column 
shown in table (1) in the variance covariance matrix. This is substituted to 
equation (26) to get p si . Equation (25) has a function to interpolate p si and p s i+1 
at two hues. The actual meaning of correction-coefficient 't* exists in p s i and p s 
i+1 . Table 2 lists p si calculated with the lightness standard deviation shown in 
Table 1 and equation (26). As shown by equation (25), p Sti has six values, and 
these six values are defined on divided hues. Other than the divided hues are 
interpolated by using interior division ratio of hue H s of an input color vs. two 
adjacent divided hues and H i+1 . For the divided hues, primary colors, i.e. 
red, green and blue as RGB digital count (d r , d g , d b ) are used, and the secondary 
colors, i.e. cyan, magenta and yellow are also used. In other words, as shown in 
Table 3, red is expressed as a hue of CRT display emitted by (d r , d & , d b ) = (255, 0, 
0). In the same manner, yellow has a hue of the CRT display emitted by (d r , d g , 
d b ) = (255, 255, 0). Relations between eight color categories except no-chroma 
colors and six divided hues are listed in Table 3. 
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TABLE 2 



Category 


ps.b 


2 


Red 


0.381 


3 


Brown 


0.912 


4 


Pink 


0.649 


5 


Orange 


0.709 


6 


Yellow 


0 


7 


Green 


0.722 


8 


Blue 


1 


9 


Purple 


0.722 


TABLE 3 



Divided hue 


Digital count 
(dr, dg, db) 


Category 


ps,b 


36.1 


(255, 0, 0) 


2 Red 


0.381 


98.2 


(255, 255, 0) 


6 Yellow 


0.000 


143.4 


(0, 255, 0) 


7 Green 


0.722 


203.2 


(0, 255, 255) 


8 Blue 


1.000 


296.5 


(0, 0, 255) 


9 Purple 


0.772 


330.5 


(255, 0, 255) 


10 Pink 


0.649 



Fig. 20 shows resultant mapping of the CRT display's cusps by giving the 
data in Table 3 to equation 19. 

In Fig. 20, line 2801 represents the CRT display's cusps, line 2802 
10 represents printer's cusps, and line 2803 does the CRT display's cusps mapped 
according to equation 19. 

In the region from red to yellow, i.e. from hue angle 0 (zero) degree to 
around 100 degree, cusp 2803 mapped agrees with cusp 2802 of the printer. 

In the range from ca. 100 degree to ca. 200 degree, i.e. the range where 
15 yellow changes to blue, cusp 2803 mapped approaches to CRT display's cusp 
2801 from printer's cusp 2802. 
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In this second embodiment, a linear mapping along the lightness axis and 
expressed by equation (18) is used for the description; however, non-linear 
mapping is applicable instead of equation (18). 

For instance, try to find the lightness — perceived by the viewer on the 
border between black and gray as well as on the border between gray and white 
also perceived by the viewer — on both the CRT display and the printer with 
respect to a color on the lightness axis (a color of L = 0). In the case of lightness 
mapping on the lightness axis from the CRT display to the printer, the following 
four points are used as reference points: (a) border between black and gray, (b) 
border, between gray and white, (c) min. lightness point on the CRT display and 
the counterpart mapped on the printer, (d) max. lightness point on the CRT 
display and the counterpart mapped on the printer. Based on these four 
reference points, overall lightness can be covered by a line graph or a curve, so 
that non-linear mapping is realized. 

Besides the mapping on the lightness axis, the present invention allows 
the colors on the CRT display's lightness axis to be mapped onto the printer's 
lightness axis. 

The categorical basic color names shown in Table 1 are used in the color 
naming by the viewer; however, other color names such as memorized colors 
which are key to determine the color reproducibility, color names set by the 
viewer per se, can be used. 

When the names of colors perceived by the viewer are used for color 
naming, the mapping can be controlled in an image gamut oriented manner, so 
that efficient and appropriate color-information-exchange is realized comparing 
with a case where a device's gamut is entirely covered. 

In this second embodiment, the variance covariance matrix in equation 
(27) is calculated through grouping the colorimetric values of the test colors 
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showed in color naming responsive to the color names given by the viewers. 
However, the variance covariance matrix can be defined by the device-drive- 
level, e.g. CRT display's drive level (d r , d g , d b ), instead of the method discussed 
above. 

Respective processes of the color information exchanging, used in the 
second embodiment, shown in Fig. 14 are detailed hereinafter. 

First, the construction of pre-mapping section 2109 is described. Fig. 21 
illustrates the structure of pre-mapping section 2109. 

In Fig. 21, source device colorimetric value input section 2901 receives an 
output from polar coordinate converter 2104 shown in Fig. 14, i.e. the lightness, 
chroma, and hue [L*, C*, HJ S of the color of the CRT display (source device). 

Source device chroma saturation calculator 2902 receives chroma C s 
supplied from input section 2901 and calculates a saturation degree of source 
device's chroma. Pre-mapped chroma calculator 2903 receives the saturation 
degree, and compresses the chroma of the color on the CRT display, then outputs 
the resultant compressed chroma to pre-mapped colorimetric output section 
2904. 

Lightness component variance database 2905 has stored the variance of 
the lightness component found to respective color names from the colorimetric 
values of the CRT display's color names. 

Cusp database 2906 for the source device has stored the lightness, 
chroma and hue of the cusp of the CRT display's gamut. 

Cusp database 2907 for the output device has stored the lightness, 
chroma and hue of the cusp of the gamut of the printer (output device). 

Cusp lightness-difference correction coefficient calculator 2908 calculates 
a cusp lightness-difference correction coefficient by using hue H s of the color by 
the color display and variance information supplied from database 2905. Cusp 
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mapping section 2909 for the source device maps the cusps of the CRT display 
onto the gamut surface of the printer based on cusp lightness-difference 
correction coefficient supplied from calculator 2908, lightness information of the 
cusps supplied from database 2906 and another lightness information of the 
5 cusps supplied from database 2907. 

An output from source device cusp mapping section 2909 is fed into gamut 
boundary lightness compressor 2910, which compresses the lightness of a color 
from the source device referring to the lightness information supplied from 
0 database 2906. On the other hand, the lightness of the color by the CRT 

^ 10 display is fed into non-chroma color lightness compressor 2911 and compressed. 
'2 ^ n output from gamut boundary lightness compressor 2910 and an 

: P output from non-chroma lightness compressor 2911 and an output from source 

m device chroma saturation calculator 2902 are supplied to pre-mapped lightness 

calculator 2912, where these outputs are converted to the lightness of pre- 
15 mapped colorimetric values. Then the converted values are supplied to pre- 
; J; mapped colorimetric output section 2904. The hue of the CRT display's color is 

O directly supplied to output section 2904. 

Next, a structure of main-mapping section 2110 shown in Fig. 14 is 
described. Fig. 22 illustrates the structure of main-mapping section 2110. 
20 Pre-mapped colorimetric value input section 21001 receives an output [L* C* 
H] u from pre-mapping section 2109. 

Gravity control point mapping section 21002 provides an output [L* C* 
H] u from pre-mapped colorimetric input section 21001 with gravity-control-point 
mapping, and supplies the resultant output to surface -control-point mapping 
25 section 21003. 

Surface-control-point mapping section 21003 receives gamut shape 
control coefficient w c supplied from gamut-shape-control-coefficient input 
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section 21004, and provides an output from gravity control point mapping 
section 21002 with surface control point mapping. The resultant output is 
supplied to the printer via output-device-colorimetric-value output section 
21005, 

Main-mapping section 2110 designs a mapping by utilizing control points 
set in the color space of the CRT display (source device) and the printer (output 
device). 

Control points are grouped into gravity control points and surface control 
points. The main-mapping includes a gravity control mapping which is 
practiced between the gravity control points, and a surface control point 
mapping which is practiced between the surface control points. All the subject 
control points on the source device and their counterparts on the output device 
are selected from the identical color category, whereby the mapping, in which 
the color category is preserved, can be realized. 

A structure of gravity control point mapping section 21002 shown in Fig. 
22 is described. Fig. 23 shows the structure of gravity control point mapping 
section 21002. Color naming database 21101 for the source device has the 
same functions as that of the first embodiment, namely, it has stored 
colorimetric values of the colors of the CRT display, and the gravity vectors for 
each color as well as variance covariance matrixes also for each color. 

Color naming database 21102 for the output device has the same 
functions as that of the first embodiment, namely, it has stored colorimetric 
values of the colors of the printer, and the gravity vectors for each color as well 
as variance covariance matrixes also for each color. 

Surface control point database 21103 for the source device has stored 
lightness, chroma, hue of the mapping control points provided to the gamut 
surface of the CRT display, or near to the surface thereof. 
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Colorimetric input section 21104 outputs the lightness, chroma, hue [L* 
C* H] u of the pre-mapped colorimetric value to the following three systems: 

The first system is to calculate a gravity difference weighing coefficient, 
and the lightness, chroma, hue [L* C* H] u are supplied to normalized distance 
calculator 21105 for the source device. 

Calculator 21105 calculates the normalized distances between the pre- 
mapped colormetric values and nine color categories based on the gravity 
vectors and variance covariance matrixes supplied from database 21101. 

Gravity mapping executing section 21106 executes mapping between the 
gravities based on the gravity vectors and variance covariance matrixes 
supplied from database 21101 as well as those supplied from database 21102. 

Gravity difference weighing coefficient calculator 21107 calculates the 
coefficient using an output from executing section 21106. 

The second system is to calculate a surface control point mapping 
restraint coefficient, and the lightness, chroma, hue [L* C* H] u are fed to surface 
control point mapping restraint coefficient calculator 21108. 

Calculator 21108 receives pre-mapped colorimetric values, the 
information of surface control point supplied from database 21103 and the 
output of database21101, then calculates the restraint coefficient. 

The third system is to use mapping executing section 21109. Gravity 
control point difference matrix calculator 21110 calculates a gravity control 
point difference matrix based on an output from pre-mapping section 21111 and 
an output from database 21102. 

Calculator 21110 receives information from pre-mapping section 21111, 
where the information contains hue, chroma, lightness of the pre-mapped 
gravity vector in the color space of the CRT display. The process by pre- 
mapping section 21111 is the same as pre-mapping section 2109 shown in Fig. 14. 
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Detailed description thereof is thus omitted here. 

Calculator 2110, on the other hand, receives information from database 
21102, where the information contains gravity vectors in the color space of the 
output device. 

Mapping executing section 21109 receives the gravity difference weighing 
coefficient from calculator 21107, also receives the surface control point mapping 
restraint coefficient from calculator 21108, and further receives the gravity 
control points difference matrix from calculator 21110, thereby executing the 
gravity-control-point mapping to the lightness, chroma, hue [L* C* H] u of the 
pre-mapped colorimetric value. 

A structure of surface-control-point-mapping section 21003 shown in Fig. 
22 is detailed hereinafter. Fig. 24 illustrates an inner structure of mapping 
section 21003. 

Surface control point database 21201 for the source device has stored the 
lightness, chroma, and hue of the mapping-control-points provided on the gamut 
surface of the CRT display or near to the surface thereof. 

Colorimetric value input section 21203 receives an output from gravity- 
control-point mapping section 21002 shown in Fig. 22. 

Source device chroma saturation calculator 21204 calculates a degree of 
source -device-chroma-saturation using the chroma component of gravity control 
point mapped colorimetric value. 

Surface control point mapping restraint coefficient calculator 21205 has 
the same function as calculator 21108 shown in Fig. 23, and calculates a surface- 
control-point-mapping-restraint coefficient. 

Surface control point difference matrix calculator 21206 calculates a 
surface control difference matrix using CRT display ! s surface control points pre- 
mapped by pre-mapping section 21207 and printer's surface control points. 
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The processes by pre-mapping section 21207 are the same as pre-mapping 
section 2109 shown in Fig. 14, thus detailed description is omitted here. 

Surface difference weighing coefficient calculator 21208 calculates a 
surface difference weighing coefficient using colorimetric values of the surface 
5 control points supplied from surface control point database 21201 for the source 
device. 

Mapping executing section 21209 receives the chroma saturation degree 
of the source device from calculator 21204, and also receives the surface control 
point mapping restraint coefficient from calculator 21205, and receives the 

10 surface control points difference matrix from calculator 21206, and further 
receives the surface difference weighing coefficient from calculator 21208, and 
further receives the gamut-shape-cntrol-coefficient from calculator 2108 in 
Fig. 14, thereby executing the surface-control-point mapping to the lightness, 
chroma, hue of the gravity-control-point mapped colorimetric value. 

15 An output from mapping executing section 21209 is supplied to the 

printer via colorimetric values output section 21210. 

Next, the processes of gravity-control-point-mapping shown in Fig. 23 and 
surface-control-point-mapping shown in Fig. 24 are described hereinafter. 

Main-mapping section 2110 converts the pre-mapped colorimetric value 

20 [L* C* H] u into main-mapped colorimetric value [L* C* H] d , i.e. final mapping 
point, by the mapping expressed in the second and third terms of equation (29). 
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where [L* C* H] u = the pre-mapped colorimetric value given by equation (16), 
matrix M = gravity control point difference matrix, 
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matrix P = surface control point difference matrix, 
W m = gravity difference weighing coefficient, 
W r = surface difference weighing coefficient, 
F = surface control point mapping restraint coefficient, 
C r = source device chroma saturation degree, 
matrix J = gamut shape control coefficient matrix. 
As expressed in equation (29), a main-mapping process by main-mapping 
section 2110 shown in Fig. 14 includes two mapping processes, namely, FMW m 
and JC r (l-F)PW r . 

Among FMW m , the MW m functions as mapping the gravity control points 
provided in the color spaces of both the CRT display and printer, and is referred 
to as "gravity control point mapping". 

Among JC r (l — F)PW r , the JC r PW r functions as mapping the surface 
control points provided in the color spaces of both the source and output devices, 
and is referred to as "surface control point mapping". 

Surface control point mapping restraint coefficient 'F' functions as keep 
balance between the gravity control point mapping MW m and surface control 
point mapping JC r PW r . Restraint coefficient 'F' interpolates the two mapping 
functions into linear model based on an input position of pre-mapped 
colorimetric value [L* C* H] u . 

The second term FMW m and the third term JC r (l-F)PW r of equation (29) 
are detailed hereinafter. 

The elements constituting the second term, i.e. gravity control point 
difference matrix M, gravity difference weighing coefficient W m , surface control 
mapping restraint coefficient F, are detailed hereinafter in this order. 

Matrix M is given by equation (30). 
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(30) 



where L* M d i? C*^ d it H% f dt i( are gravity vectors of color category f i' in the 
printer. 

The color category employs the nine color names shown in Table 1. The 
gravity vectors are found in this way: Sample data in color naming are 
converted into values a j3 r with equation (28), then gravities of respective color 
categories are found by the converted values in equation (31), and they undergo 
polar coordinate conversion. 




In equation (30), L* Mtd>i , C% d i? H% di , represent pre-mapped gravity vectors, 
where the gravity vector of color category T of the source device is mapped with 
equation (16). The color category employs the nine color names listed in Table 
1, and the gravity vectors are found by extracting samples 1 gravities from 
equation (31), then they undergo polar coordinate conversion, this is the same 
procedure as that in the case of the output device. As equation (30) shows, 
when the pre-mapped colorimetric values [L* C* H] u are extracted from gravity 
vectors of the source device, matrix M is to map the pre-mapped colorimetric 
values [L* C* H] u to gravity vector [L* C* H]^ of the output device. Therefore, 
matrix M has a function of mapping the gravity vector in the color space of the 
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source device to that in the color space of the output device. 

Next, gravity difference weighing coefficient W m is described hereinafter. 

Weighing coefficient W m weighs mapping control over matrix M with 
respect to each gravity of respective colors based on an input position of the pre- 
mapped colorimetric [L* C* H] u , thereby realizing the mapping control over 
matrix M with respect to any pre-mapped colorimetric value [L* C* H] u . 
Coefficient W m is given by equation (32). 
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where D ti represents an intermediate normalized distance. Coefficient W m 
increases at smaller intermediate normalized distance D t i( and decreases at 
greater intermediate normalized distance D t ; , which is given by equation (33). 

Dt =VDs 
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where D s represents a color categorical normalized distance vector of the source 
device's colorimetric value. Gravity mapping scaling matrix f V' obtains 
25 intermediate normalized distance D t by mapping color categorical normalized 
distance D s of the source device's colorimetric value to the color space of the 
output device. 
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On the other hand, color categorical normalized distance D i is the 
distance from a certain colorimetric vector X = [ a )3 7 ] t to the gravity of color 
category T and given by equation (34). 

Di = *s](X - (34) 

where 11 i represents the gravity vector of color category 'i 1 and E 4 represents 
the variance covariance matrix of color category 'i\ 

As shown in Fig. 19, a spread of color category's distribution depends on 
respective color categories, and a distance in CIELAB space is not proportionate 
to a color categorical distance. In other words, red category (A) has a smaller 
spread than green category (A), and a distance in CIELAB space seems a 
greater distance for the red category from the color categorical view of point. 
Color categorical normalized distance D t is a normalized distance with 
consideration of distribution-spreads of respective color categories, and distance 
Di is proportionate to a color categorical perspective measure. 

Gravity mapping scaling matrix f V f in equation (33) is described 
hereinafter. This matrix 'V is given by equation (35). 
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where $ represents a gravity mapping scaling matrix, and this matrix has a 
function of converting a distance-relation between gravity vectors in the CRT 
display into a distance-relation between gravity vectors in the printer. 
Geometrically, it maps the gravity vectors of the source device to the gravity 
vectors of the output device. 

Gravity mapping scaling coefficient matrix Ohas element v lm . This 
element v lm is a scaling coefficient which converts normalized distance D s lm — 
extending over the gravity of color category T to the gravity of color category 'm' 
in the color space of the source device — into normalized distance D d ^ — 
extending over the gravity of color category T to the gravity of color category 'm' 
in the color space of the output device. 

Matrix <E> has all the scaling coefficients with respect to every possible 
combination of nine color categories, thus the gravity vectors in the color space 
of the source device are mapped to the gravity vectors in the same color category. 

However, as shown by equation (35), scaling coefficient (v n , , v mm ) for its 

own distance becomes zero (0). 

Geometrically, matrix $ means to map the nine gravity vectors of the 
source device to those of the output device. From the viewpoint of color 
reproducibility, it corresponds to the mapping between the original and its 
counterpart having the same color category by using the normalized distance. 

Intermediate normalized distance D t i is a distance found by mapping the 
color categorical normalized distance in the source device to the output device 
with matrix 'V*. This mapping corresponds to the mapping of the gravity of the 
same color category. 

Next, the way of finding weighing coefficient w v i is described with respect 
to gravity-control-point-difference matrix 'IVT of equation (35) by using the 
normalized distance mapped to the output device's color space. Weighing 
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coefficient w ^ is given by equation (36). 
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5 (36) 

if D s, i = 0 then w v? i = 1, w v k ^ f = 0 

if £>s,k* i - 0 f/ze/z w Vf jt „ f = 1, w Vj , ^ f . „ * = 0 

Equation (36) functions as covering the mapping of gravity to the places 
10 other than the gravity, and interpolates gravity's mapping information 
responsive to the positional relation between color vector X and gravity vector 
U before carrying out the function. 

In equation (35), column vector element [v^, v 2i , , v 9 J of gravity 

mapping scaling coefficient matrix <E> is a scaling coefficient in color category Y, 
15 and these scaling coefficients are weighed by weighing coefficient w^ . 

Weighing coefficient w vi is defined by B/A, where 'A 1 represents the sum of 
inverses of the normalized distances between the gravities of nine color 
categories, and 'B f represents the inverse of the normalized distances between 
the gravities of color category ! i\ Therefore, the gravity of the color category 
20 nearer to color vector X has more important scaling coefficient, and the gravity 
of the color category farther from color vector X has less important scaling 
coefficient. 

In summary, pre-mapped colorimetric value [L* C* H] u is converted to 
color categorical normalized distance D s with equation (34), then converted to 
25 intermediate normalized distance D t with gravity-mapping-scaling-matrix 'V of 
equation (33). Then in the output device's color space, gravity-difference- 
weighing-coefficient w m is calculated with equation (32). 
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Surface control point mapping restraint coefficient £ F is described 
hereinafter. As mentioned previously, coefficient T' functions as keeping 
balance between gravity-control-mapping MW m and surface-control-point- 
mapping JC r PW r , and based on pre-mapped colorimetric value [L* C* H] u , 
5 coefficient T* interpolates these two mapping functions in a linear manner. 

In other words, when restraint coefficient T' takes T (one) in equation 
(29), the third term, i.e. surface control point mapping JC r PW rJ does not work at 
all, and only the second term, i.e. gravity control point mapping MW m , does work. 
On the contrary, when coefficient ( F takes c 0' (zero), only mapping JC r PW r works 
10 and mapping MW m does not work at all. 

The surface control point is used in order to set restraint coefficient T\ 
The surface control point is disposed on the device's gamut surface or near to the 
surface. The surface -control-point functions as a control point for mapping the 
gamut boundary. 

15 In this second embodiment, eight surface-control-points are defined for 

the simple description; however, a number of the control points can be selected 
arbitrarily. 

The eight surface control points are defined in digital counts of the device 
as follows: the primary color: red R:(d r , d g , d b ) = (255, 0, 0), green G: (d rJ d g , d b ) = 
20 (0, 255, 0), and blue B: (d r , d g , d b ) = (0, 0, 255). Secondary colors: yellow Y: (d r , d g> 
d b ) = (255, 255, 0), magenta M: (d r , d g , d b ) = (255, 0, 255), cyan C: (d r , d g , d b ) = (0, 
255, 255), white W: (d r , d g , d b ) = (255, 255, 255), black B: (d r , d g , d b ) = (0, 0, 0). 
Restraint coefficient 'F is given by equation (37). 

25 F=w u f 1 + w £2 f 2 + • • • + Wf9 f 9 (37) 



where f { and w f ;i are given by equations (38) and (39) respectively. 
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'/ As,/,;=0 then w f • = 1, w f tk „ , « 0 

*/ ^>s, /,* „ / = 0 tfze« W/ k „ l = 0, w ffl „ ; ^ k = 0 

15 E sfi represents an Euclidean distance between gravity vector Ui of color 

category Y and color vector X in the CIELAB space. E sai represents an 
Euclidean distance between surface control point T and color vector X in the 
CIELAB space. Both of E sfi and E sa i show distances in the CRT display's color 
space. 

20 When E s f i = 0, namely, color vector X agrees with gravity vector tt it fi = 1 

which means constraint force becomes maximum extent, then surface -control- 
point-mapping JC r PW r does not work at all, and only gravity-contorl-point- 
mapping MW m works. On the contrary, when E sai = 0, namely, a test color 
agrees with the surface control point, f | = 0 which means no constraint force is 

25 available, then gravity-control-point-mapping MW m does not work at all and 
only surface-control-point-mapping JC r PW r works. 

In equation (39), D sfi represents the normalized distance between gravity 
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vector Ui of color category T and color vector X in the source device's color 
space. 

When color vector X agrees with gravity vector H x of color category T, w f i 
= 1 (one), and weighing coefficient w fk!=1 for other than color category V becomes 
5 0 (zero). Accordingly, w fi is a weighing coefficient for constraint force £ 
determined following the position of color vector X. Nine gravity vectors have 
constraint forces f, respectively. Those forces are weighed by w ti and then 
reflected to restraint coefficient T\ 
P The third term of equation (29) has four elements, i.e. surface control 

h >4 10 point difference matrix T\ surface weighing coefficient w r , source device's 
% g chroma saturation C r , and gamut shape control coefficient matrix 'J\ and these 

fTj four are described hereinafter in this sequence. 

=31 

'~ Gravity control points of source device's color space are mapped to output 

device's color space by gravity control point mapping FMW m . In the same 
15 manner, the surface control points in CRT display's color space are mapped to 
those in printer's color space by surface control point mapping JC r (l - F)PW r . 
In other words, elements of matrix 'M' replaced with the surface control points 
are matrix T', and matrix T' is given by equation (40). 
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(40) 



When surface control point mapping restraint coefficient T' .= 1, 
coefficient 'F restraints the entire mapping between the surface control points. 
25 Surface difference weighing coefficient w r is calculated with equation (41) based 
on the positional relation between pre-mapped colorimetric value [L C H] u , and 
weighs the mapping operation by matrix ( P\ 
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(41) 



Saturation C r is given by equation (17) and functions as a chroma- 
dependent weighing coefficient such that J(l - F)PW r mapping works mainly in 
a high chroma region. Therefore, the third term JC t (l - F)PW r does not work 
at all on the lightness axis. 

Matrix 4 J' functions as mapping the gamut boundary of the source device 
to that of the output device with equation (42). 



I=£ 15 
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1 0 0 
0 w c 0 
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(42) 
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where gamut shape control coefficient w c is given by equation (43), and it only 
controls chroma 'C*' and maps CRT display's gamut boundary to printer's 
gamut boundary. 



_ d, max, t ~ ^ c,s, max,f ~ FM c w m 

C r (l-F)P c w r 



(43) 



where C* d max t corresponds to black point 21301 in Fig. 25, and black point 21302 
25 is a point [L* C*H] CS _ maxd mapped with equation (44). Fig. 25 illustrates chroma 



c,s,max, 



d and C* d max t when coefficient w c is found. 
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Equation (44) re-maps the point [L* C* H] csmaxu which has been pre- 
mapped from the point [L* C* H] csmaX(d with equation (16). The point [L C 
H] C(Stmax;d is main-mapped with equation (29). When black point 21302 moves to 
point 21301, a mapping point is determined to be on the gamut boundary of the 
output device. Therefore, equation (45) is established and control coefficient w c 
is found therefrom, then an adjustment coefficient, which maps the gamut 
boundary of the source device to that of the output device, can be figured out. 



15 



20 



25 



C* c , s , max>t = C* c , s>maX(U +FM c +w c C r (l-F)PW r 



(45) 



In equation (45), vector M c is a one-row, nine-column type row-vector, which is 
formed by taking out a chroma component out of gravity control point difference 
matrix 'M\ Vector P c is a one-line, eight-column type row-vector, which is 
formed by taking out a chroma component out of surface control point difference 
matrix T\ 

As discussed above, gamut shape control coefficient w c functions 
compressing the chroma component of C r (l - F)PW r to a linear type so that a 
color on source device's gamut boundary can be mapped onto the gamut 
boundary of the output device. 

Exemplary Embodiment 3 

A color information exchanging device in accordance with the third 
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exemplary embodiment of the present invention is demonstrated hereinafter 
with reference to Fig. 26. 

Fig. 26 is a block diagram of a color management system introducing a 
color information exchanging device in accordance with the third exemplary 
embodiment of the present invention. 

This third embodiment is designed based on an assumption that the 
gamut of a source device is not so much different from that of an output device. 
Therefore, a case, where the source device is a CRT display and the output 
device is a projection display, corresponds to the third embodiment. Another 
case, where both of the source and output devices employ CRT displays and yet 
the color temperatures of respective reference white color points are different, 
also corresponds to this third embodiment. 

In the following description, the source device is a CRT display having a 
color temperature of 9300K (hereinafter referred to as 9300K-CRT display), and 
the output device is also a CRT display having a color temperature of 5000K 
(hereinafter referred to as 5000K-CRT display). It is assumed that an 
appearance of an image reproduced on the 5000K-CRT display is to agree with 
the image on the 9300K-CRT display, where a white color point appears more 
bluish. 

In Fig. 26, source device driving unit 21402 displays an image on the 
source device 21401 (9300K-CRT display), and also sends out colorimetric values 
of displayed colors on the CRT display to color-information-exchanging-device 
21400 by referring to source device's input-output characteristics stored in 
profile 21403 for the source device. 

Output device driving unit 21410 drives output device 21412 (5000K-CRT 
display), and also converts the colorimetric value supplied from exchanging- 
device 21400 to a drive signal which can be appropriately displayed on the 
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5000K-CRT display by referring to input-output characteristics of the 5000K- 
CRT display. The input-output characteristics are stored in profile 21411 for 
the output device. 

Color information exchanging device 21400 of the present invention is 
5 described hereinafter. Polar coordinate converter 21404 converts the 
colormetric values of the displayed color on 9300K-CRT display into lightness, 
chroma, and hue. These converted values are supplied to two systems. 

One of the systems is a gamut mapping system. In this system, after the 
process in mapping processing section 21408, gamut mapping is practiced to the 
10 lightness, chroma and hue of 9300K-CRT display. 

The other system is a system for finding a gamut shape control coefficient. 
In this system, the lightness, chroma and hue of 9300K-CRT display undergo 
the processes by max. chroma color setting section 21405 and max. chroma color 
mapping section 21406, then the lightness, chroma and hue are supplied to 
15 gamut shape control coefficient calculator 21407, where an adjustment 
coefficient is calculated as a gamut shape control coefficient which can map the 
gamut boundary of 9300K-CRT to that of 5000K-CRT display. 

An output from mapping section 21408 is converted to a colorimetric 
value of orthogonal coordinate from the lightness, chroma and hue of the polar 
20 coordinate by orthogonal coordinate converter 21409, then the converted value 
is supplied to output device driving unit 21410. 

In the mapping from 9300K-CRT display to 5000K-CRT display, the color 
information exchanging device of the present invention advantageously makes 
the color appearance of the image displayed on 5000K-CRT display agree with 
25 that on 9300K-CRT display by determining mapping points so that the color 
category of the source device agrees with that of the output device. 

This third embodiment differs from the second embodiment in no pre- 
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mapping process. Other than this difference, the structure and process are the 
same as those of the second embodiment, thus the descriptions thereof are 
omitted here. 

In this third embodiment, since both the source and output devices 
5 employ CRT displays, gamut mapping due to the difference of the gamut shapes 
is not needed. The mapping in this third embodiment corresponds to the 
mapping, which produces a pair of which color appearances agree with each 
other, between the gamuts of almost the same shapes. 

□ Therefore, the pre-mapping process, which corrects the difference in the 
r] 10 gamut shapes, is not needed, and the mapping of the third embodiment makes 
;~ color appearances agree with each other for realizing color reproduction in a 

natural way. 

y " In conclusion, the colormetric values of the test colors or device drive 

j"* levels obtained by color naming done by at least one viewer are grouped into 

5 

U 15 color categories, where the color naming is performed to a plurality of test colors 

□ rendered on the source and output devices. Variance and other statistical 
characteristics are used for defining a normalized distance of each color category, 
and the relation between color perception of the viewer and color reproducibility 
is understood quantatively, so that gamut mapping is designed. This is an 

20 advantage of the present invention. 

The present invention produces another advantageous effect: The color 
spaces CIELAB of the source and output devices are grouped into a plurality of 
color categories. Each device produces a normalized color space that is 
normalized by using a gravity vector calculated for each color category of 

25 respective devices and a variance covariance matrix expressing extension from a 
gravity point as a center. Then control points of the same color categories 
provided in the normalized color spaces are mapped, whereby color information 
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is exchanged between the two normalized color spaces by agreeing the color 
categories with each other and preserving color perceptional quantity including 
chroma of two colors in respective normalized color spaces. This arrangement 
allows the entire impression of the images rendered on respective devices to 
appear as little difference as possible. 

There is still another advantageous effect of the present invention: Apre- 
mapped color is produced by mapping the color rendered on the source device, 
where the lightness and chroma of this color have been changed. Then the 
lightness, chroma and hue of the pre-mapped color are changed by using a 
weighing coefficient based on a normalized distance for re-mapping the pre- 
mapped color. This arrangement allows the mapping between the devices 
having different gamuts to be designed as a gamut mapping based on clear 
references found through a standpoint of agreement of color categories. 

When the gamut shapes of the source device and the output device are not 
so much different, the pre-mapping process, which maps source device's gamut 
into output device's gamut, can be omitted. Then a mapping which realized an 
agreement of color categories is performed, so that the color appearance of the 
source device agrees with that of the output device. In other words, this 
arrangement can advantageously realize a matching of color appearances. 



